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Abstract 

With the aid of a Fermi- Walker chart associated with an orthonor- 
mal frame attached to a time-like curve in spacetime, a discussion is 
given of relativistic balance laws that may be used to construct models 
of massive particles with spin, electric charge and a magnetic moment, 
interacting with background electromagnetic fields and gravitation de- 
scribed by non-Riemannian geometries. A natural generalisation to 
relativistic Cosserat media is immediate. 

1 Introduction 

The language of relativistic fields on a spacetime offers, probably, the most 
succinct and complete description of the basic laws of physics. From these 
laws we expect to accommodate the enormous amount of data derived from 
experiment, predict new phenomena and modify or reject laws in accordance 
with observation. Of necessity this process invariably involves some level of 
approximation. As experimental technology improves it becomes possible 
to refine such approximations and thereby seek more detailed experimental 
verification of theoretical predictions. 

A field description of relativistic gravitation was formulated by Einstein 
almost a century ago. Since then it has proved one of the most successful 
descriptions of large scale gravitational phenomena in all of physics and has 
profound implications for the notions of space and time on all scales. Nev- 
ertheless this theory of gravitation cannot be complete and it poses many 
challenges in both classical and quantum domains. Modern experimental 
tools of increasing sophistication are now being brought to bear on some of 
these problems with the search for classical gravitational waves, the break- 
down of the "equivalence principle" and the influence of gravitation on the 
"collapse of the wave- function" . Along with these technological advances 
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new theoretical avenues continue to be explored in which all the basic inter- 
actions between matter are subsumed into a coherent framework. 

Viable approximation schemes that link the results of experimental probes 
and theoretical descriptions of the basic interactions will remain as bridges 
between these continuing developments. Since the early days of relativity 
it was recognised that the derivation of the motion of extended "matter" in 
gravitational fields was a non-trivial problem and great effort was expended 
in deriving effective approximations. The work of Einstein, Infeld, Hoffman, 
Mathisson, Papapetrou, Pryce , Tulczyjew and others ^ HI 01 S culmi- 
nated in a series of papers by Dixon, [HI El EH EI , Ehlers and Rudolph 
jl4j and Kunzle jS], where it was demonstrated that a viable scheme could 
be established, given suitable subsidiary conditions, in which an exact field 
description could be replaced by an approximation procedure in which ex- 
tended matter was regarded at each order of approximation by a finite set of 
(non-unique) matter and electromagnetic multipoles. Such a scheme offers in 
principle a viable approach to the motion of extended matter and deserves 
to be taken seriously in the context of modern experiments that purport 
to analyse the fundamentals of the gravitational interaction with spinning 
matter. However theories of gravitation with spinning matter sit somewhat 
un-naturally in the pseudo-Riemannian (i.e. torsion free) environment in 
which the whole Dixon framework was formulated. Furthermore it is not 
necessary to have "spinning sources" to accommodate gravitational fields 
with torsion. Indeed a Lorentz gauge-covariant formulation of the Brans- 
Dicke theory naturally yields a metric-compatible connection on the bundle 
of linear frames over spacetime with torsion determined by the gradient of 
the Brans-Dicke scalar field ^5]. Since many modern variants of this theory 
(and various "string theories") abound with scalar gravitational fields one 
should seriously consider the possibility that gravitation may have a tor- 
sional component that is absent (by hypothesis) in the pseudo-Riemannian 
framework given by Einstein. 

The multipole analysis formulated by Dixon relies heavily on the use 
of bi-tensors and tensor densities and the pseudo-Riemannian structure of 
spacetime. To extend this methodology to general spacetime structures a 
new approach is adopted here based on the use of differential forms, structure 
equations, a (generalised) Fermi chart associated with a timelike worldline 
and the transport of a Fermi- Walker frame along spacelike autoparallels of 
a metric-compatible connection with torsion. This approach may be con- 
trasted with the use of adapted coordinates employed by Tulczyjew and 
Tulczyjew |13| in a torsion-free context. Although they employ a Fermi 
chart in discussing the equations of motion of an electrically neutral parti- 
cle, their methods, detailed calculations and conclusions differ significantly 
from the results to be discussed below even in the case considered, where 
the torsion is zero. In particular their resulting equations admit unnatural 
solutions in Minkowski spacetime as will be discussed in section VIII. 
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The approach to be considered here considers the modelhng of rela- 
tivistic massive particle interactions witli background electromagnetic and 
gravitational fields with possible torsion. As stressed above it is assumed 
that the underlying laws describing all phenomena can be formulated in 
terms of tensor and spinor fields over a manifold. Depending on the phys- 
ical scales involved such fields are in general subject to either the laws of 
classical or quantum field theory. In this paper emphasis is on classical 
fields. In particular it is assumed that classical gravitation can be described 
in terms of a geometry on a spacetime and the electromagnetic field is a 
closed 2-form on this manifold. Like Newtonian gravitation, electromag- 
netism is distinguished from other gauge fields by its long range character. 
It can also be related to the curvature of a connection on a principle U{1) 
bundle over spacetime. The connection on the bundle of linear frames over 
spacetime will be that induced from the bundle or orthonormal frames. 
Thus the local gauge group for gravitation may be regarded as the covering 
of the Lorentz group 50(1,3). Shorter range 2-form (internal symmetry 
non-Abelian Yang-Mills) fields, regarded as curvatures associated with con- 
nections on other principle group bundles are also responsible for the basic 
interactions in Nature. Matter fields (regarded as sections of associated 
bundles) carry representations of these Yang-Mills groups and their gauge 
covariant interactions with gauge fields are responsible for the "charges" of 
massive field quanta. All fields influence the geometry of spacetime and in 
turn the gravitational field influences the behavior of all other fields. The 
field laws for the coupled system of gauge fields and matter may be com- 
pactly derived as differential equations that render extremal the integral of 
some (action) 4-form on spacetime. This integral is constructed to be in- 
variant under the action of all local gauge groups. The resulting classical 
variational field equations are thereby covariant under changes of section 
induced by associated local changes of frame (gauge). In this article the 
Yang-Mills internal symmetry group made explicit will be restricted to elec- 
tromagnetism. It is the gauge invariance of some action integral that is of 
paramount importance in the developments below. 



2 Balance Laws 

Classical point particles arc an abstraction. In general they offer a useful 
approximation to the behavior of localised matter and form the basis of rigid 
body dynamics in Newtonian physics following the laws of Newtonian parti- 
cle dynamics. The behavior of Newtonian deformable continua require new 
laws that exploit the high degree of Killing symmetry in spacetime devoid 
of gravitation where matter, (observed from weakly accelerating reference 
frames) has a small speed relative to that of light. The balance laws of en- 
ergy, linear and angular momentum must be separately postulated in order 
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to provide a framework for the description of non-relativistic deformable 
continua. However the evolution of such continua can only be determined 
from such balance laws once a choice of constitutive relations (Cauchy stress 
tensor) consistent with the laws of thermodynamics has been made. Such 
relations serve to render the balance laws predictive and are either phe- 
nomenological in nature or arise as coarse-grained approximations applied 
to a statistical or quantum analysis that recognises the atomistic nature of 
matter. Newtonian balance laws refer to arbitrary distributions of matter in 
space. Point particle and rigid-body dynamics arise as zeroth and first or- 
der multipole approximations jl^l- There are other types of approximation 
that are useful when the continuum is slender in either one or two dimen- 
sions relative to the third (i.e. rods and shells). Furthermore the internal 
state of such structures can be expressed as degrees of freedom associated 
with a field of frames in the medium (Cosserat models). In all such cases a 
successful phenomenological description is sought in terms of a small num- 
ber of physical parameters (mass, moment of inertia tensor, elastic moduli, 
pressure, density, thermal and viscosity coefficients...) that is amenable to 
experimental verification. For extended structures, such as fluids or localised 
structures such as strings or rods it is the relatively weak interaction of the 
continuum with its environment that renders such a description effective 
and enables one to distinguish externally prescribed forces and constraints 
from the responses that they produce. 

An analogous strategy is sought to describe the relativistic motion of 
localised matter from an underlying field theory involving gravitation and 
electromagnetism. Since the matter fields may carry representations of the 
Lorentz group (or its covering) as well as U{1) one expects to accommodate 
electrically charged particles with some kind of spin. The basic approxima- 
tion considered here treats the gauge fields and gravitation as a background 
that is prescribed independently of the matter fields. These on the other 
hand are required to satisfy gauge-covariant field equations in the presence 
of the prescribed gauge fields and gravitation. The metric of the back- 
ground spacetime will, however, possess no particular isometries to provide 
the analogues of the Newtonian balance laws. Nevertheless the principle of 
local gauge covariance for the interactions between fields provides powerful 
dynamical constraints on the field variables that arise in any variational de- 
scription. This observation is well known and its implications can be found 
in many parts of the literature. Such constraints will be explored below in 
the context of adapted orthonormal frame fields and associated coordinates 
based on a Fermi chart associated with a time-like curve in spacetime. 

The equations of motion proposed by Dixon offer a consistent dynam- 
ical scheme ^1] for the classical behavior of ""spinning matter" I17j 
in pseudo-Riemannian spacetimes. In the lowest pole-dipole approximation 
they can be solved analytically for electrically neutral particles in gravita- 
tional fields with high symmetry jJH] but in general recourse to numerical 
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methods is required ^HI, [201 We shall recover these equations below 
in spaces with zero torsion albeit without recourse to bi-tensors and ten- 
sor densities. Furthermore we shall show that in background spacetimes in 
which gravitation is described by geometries having connections with tor- 
sion, electrically charged spinning particles deviate from such histories due 
to prescribed interactions between the background fields and the spin and 
charge of the particle. 

We argue that in more general spacetimes even electrically neutral spin- 
less particles may have histories that follow autoparallels of connections with 
torsion instead of geodesies of the spacetime metric. 



3 Fermi Coordinates 

Let M be a smooth manifold with a metric tensor g and a metric compatible 
connection V. Denote the tangent space to M at point p by TpM. For any p 
and any vector v in TpM let 7„ : A i— > 7,;(A) € M be the unique autoparallel 
(V^„/7^,' = 0) of V with 

lv{0)=p (1) 



7'M = v. (2) 

Provided 7t,(l) exists denote it by ex.pp{v) i.e. 

expp(t>) = 7^1) G M. (3) 

The exponential map exp^ : TpM i— > M is therefore defined in a neighbour- 
hood of G TpM and is in fact a diffeomorphism |22l [22] . 

Suppose M is a four dimensional spacetime and g has Lorentzian sig- 
nature (— , -|-, -|-, -|-). The geometry of M is described by g and V which is 
metric-compatible but not assumed to be torsion free. Let a :t ^ i-^ M 
be a time-like future pointing affinely parametrised curve {g{&,a) = — 1) 
and J- = {Xq, Xi, X2, X^} a gi-orthonormal frame on a with a{t) = Xo[(j{t)- 
Then a Fermi chart ^ with Fermi coordinates (t, ) relative to the 

curve fj, frame J- and arbitrary point 17(^0) on the image of a is defined by: 

^'0(exp,(t)t;) = t-to (4) 



'^Hew.(t)v) = x'' (5) 

where v = Yl^=i^'^ -^j ^ '^a{t)M- In the following repeated Latin indices 
will be summed from 1 to 3 while repeated Greek indices will be summed 
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from to 3. The chart is defined on any open neighbourhood U of a con- 
taining points that can be joined to a by an affinely parameterised unique 
autoparahel meeting a orthogonally. Such a definition gives rise to a class 
of Fermi charts for a given curve a with frames related to each other by 
local S0{3) transformations and translations of origin a"(to) along a |24j . 
In any particular chart with coordinates {t,x^} one has the local vector 
fields {dt,dk = gfr} on U which by construction are g-orthonormal on a. 
Furthermore, by definition, the space-like autoparallels that leave a at 
{t, 0, 0, 0} can be affinely parameterised in this chart by the equations 

t-to = T (6) 



= /A (7) 

where (p^)^ + (p^)^ + (p^)^ = 1. Henceforth the origin is fixed with to = 
0. The parameters are the direction cosines of the tangent to such an 
autoparallel, with respect to T at {t, 0,0,0}. It will prove useful below to 
introduce Fermi-normal hyper-cylindrical polar coordinates {r, A,p'^} where 
{p^ p^) = 1 on part of U. These are naturally related to the coordinates 
{t, x^} by the transformation t = t and x^ = p^X. Thus any function / of 
can be regarded as a function F of t, X,p^ ,p'^ ,p^ on a suitable 

domain 

fit, x\x\x^) = fit, \p\ \p\ Ap3) = F(r, \,p\p\p^). (8) 
Using the inverse relations 

\^ix\x^,x^) = xJxJ (9) 



p\x\x\x^)= 3 (10) 

Xyx^ ,x^ ,x-^) 

where A 7^ one has the directional derivative 

dx = p'^dk (11) 



and the relations dxp^ = and = — ^ for j k. If F is analytic in A: 



00 ,f, 

FiT,\,p\p\p^) = Y,^{d\F)iT,Q,p\p\p^) (12) 



A;! 

A:=0 



E E )jP''---P'Hdn.--dU){t,^M (13) 



fc=0 ii,i2...jfc=l 
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Such expansions will be used to represent the structure 0-forms and 1-forms 
to be introduced below and partial derivatives of F with respect to A will 
be denoted F'. To eliminate excessive notation, expressions evaluated on 
the curve a will be understood to be represented in the coordinate system 
{t, x'^} and the abbreviation / often used for / restricted to a{t). 



4 Ortho-normal frames on U 

The coordinate frame = {dt ,dk} on U is only orthonormal on a where 
it coincides with T. However one can displace the frame by parallel 
transport along the space-like autoparallel curves that emanate from a using 
the connection V. Since this connection is assumed metric-compatible this 
process will give rise to a field of orthonormal frames O onU jJS] . The next 
task is to calculate the orthonormal coframe field {e°} dual to O in terms 
of the coframe {dt, dx^} dual to JF. If one assumes that the curvature and 
torsion are analytic in the vicinity lA oi a this can be done by recursively 
differentiating the structure equations with respect to the radial variable A 
along a fixed spacelike autoparallel 7^,^ on U. 

At each step of the process one evaluates the (derivatives of) structure 
forms at A = thereby building up a Taylor series in A for each. To initialise 
the process one must fix the zeroth order terms in the expansion. The 
coframe field {e°} along 7^^^ is defined to satisfy 

Vy^^e- = (14) 

i.e 



ie.oj^p = (15) 

in terms of the interior derivative of the connection 1-forms of V in the 
coframe O. Since all points of lA can be connected by such autoparallels 
we conclude that must be independent of (i A in the frame O and in 
particular uj°-p{dj) = on it. Since is dual to Xj = dj then e* = dx^. But 
dx* = dA + A dp' so e' pulls back to dA at A = 0. Also since o" = ^ 
then eP = dt. In general 

e" = e"(5t) dt + e"(aA) dA + 6^(7^) d/. (16) 



But along 7^^, dx{e'^{dx)) = VQ^{e°^{dx)) = 0, since e" is parallel and 
79^ is an autoparallel. Again, we may conclude that on U the structure 
function e°'{dx) is independent of the coordinate A. It follows that e'^{dx) = 
e^{dx)\a = and e^{dx) = e^{dx)\a = dt{dx)\a = 0. Thus one may write: 



:^ = dt + U^ dX + £^ (17) 
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where J^^ are forms dependent on t, A, and 

no = 0, = / (18) 



s'^ = 5^y + 5'ji or (19) 

with 1-forms 3^, Q*^ independent of di, d \ but components dependent on 
t, A, p'^ . Thus: 

e° = dt + yjdp> (20) 



e'' = J^'' dt + p'^dX + Q'' jdjp. (21) 
where 

= 1 (22) 

= 2*= = jr'^ = 0. (23) 

Similarly, since oj°'j3{dx) = on W, oj°'p{dj) = 0. Using ^ = and 

dp> ^ dx^ p^ dx^ ^ ^ 

one may write 

= ^"^ dt + C"^ (25) 

where A^p are forms dependent on t, A, p'^ and C^p are 1-forms inde- 
pendent of d d A but with components dependent on t, A, p^ . 
Furthermore 

c^"/3(5t)U(t)=i"/3 (26) 

= C% = u;-;3(9,.) dp^' = 0. (27) 

The structure functions can be determined in terms of properties 
of the frame P on a. Suppose the elements of P satisfy: 

F^X^ = nJ{T)Xp (28) 
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where Qai3 = —^/Sa, ^oj = and = Cij^ ^ki'T) in terms of the alternating 
symbol e^^^, indices on Q and e being raised or lowered with the Lorentzian 
metric components t] = diag{— 1,1,1,1). The (generalised) Fermi- Walker 
derivative j26j above is defined by: 



F^Xa. = V<^1« + g{&, Xa)Aa - giAa,X^)& (29) 

= 6j^'a{cr)Xf, + g{&,Xa)Aa - g{Aa,Xa)a (30) 

where Aa = VctCJ. Such a transport law maintains Xq tangent to a and 
J- g-orthonormal even though a need not be an autoparallel of V. So given 
a = Xo = dt with e''(X^) = d'j^ apphcation of e" to (EHl) yields: 

u^idt) = Qa" - ma e%A„) + g{A.,X^) 5^ (31) 

or 

(^''i{dt)=g{A^,Xi) (32) 

Cj'j{dt) = ni\T) = ni,{T). (33) 

Thus 

i"^ = Q:^p{dt) ^ (34) 

say, is fixed in terms of the orthornormal components of the acceleration 
Aa of a and the instantaneous angular velocity fifc of the frame J- along a. 
Since a is normalised and time-like the acceleration is orthogonal to a. It is 
convenient to denote the acceleration by A and write 

A,=g{Aa,X,). (35) 

If one chooses Vt^ = Q then T is said to be a non-rotating frame along o. 
Thus with = and A"^ p = A"fs the equations (11^ subject 

to (P7jl . can be substituted into the structure equations 

[271 im 123] that define the torsion and curvature of V: 

de" = Ae^ + T° (36) 
do;"^ = Aa;^/3 + ^"/3 (37) 
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where = ^ T"^i, l\e^ and i?"/3 = ^ py^v Ae^ are the torsion and 
curvature 2-forms in the frame O. Equations (|36|) . (|37|) must be satisfied for 
ah A along auto-parahels that leave C7(i) in any direction so by equating 
forms containing dXA dt and dXA dp^ on each side one can find differential 
equations for the structure forms , 3^, Q* , A°'p , C^p : 

^°' = ^°fc/ +r% (38) 

y = C\ / + T\o / 3^ + p'^ Q'^ (39) 
Q''= dp' + C\ / + r fco / 3^ + / Q'^ (40) 
A'p' = K'p^yP^T^ (41) 

0' = ^Vo/3^ + ^VnP''Q"- (42) 
These have unique solutions satisfying the conditions: 

P = 1 (43) 

= 0, q!" = 0, C'^p = (44) 

i"/3=A°/3. (45) 

By successively differentiating the differential equations above with re- 
spect to A and evaluating the results at A = one may construct the Taylor 
series representation of each analytic structure form about A = in terms 
of (derivatives of) the components of the curvature and torsion tensor and 
the transport functions A"^ on a. The first few iterations of this process 
can be readily established. To first order in A one finds: 

^"'=/(A"fc +f-fco) (46) 
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y' = (47) 

&'=dp' (48) 

i"/3' = ^Vo/ (49) 

C"^ ' = (50) 

Radial derivatives of the curvature and torsion start to appear at the 
next order: 

" = i?"feno;>>' + T"feo + T"fc^ /p" (A^ + T^o ) (51) 

y" = T^knp''dp^ (52) 

&" = f'knP'' dp^ (53) 

i"/' = ^Vo'/ +i?V^/;>"(A^ +f'^„o) (54) 

C%" = EVn/dp"- (55) 



At the third order: 

=i?"fc,0'p>' + i?"fcn^ P>'p" (A^™ + f ) + 
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3^'" = i?VnP'"p'rfp" +2fV/dp" +f\^f\,p^p^ dp'^ (57) 



"' = R\nq d + 2f\n 'p'' dp^ + T f\, p'' dp'^ (58) 

+ i?"/3fcM ^''ngO P" P^ + ^^0 'p' P" + (59) 

+ ^rpk^. f^'mu P^p'^'f' (A\ + f ^0 ) 

C"/" = 2i2Vn'/ dp"" + R'^pui f'^qp^p^ dp'^ . (60) 

In these expressions the directional derivatives d"^x for the appropriate 
order m are to be expressed in terms of derivatives with respect to x^: 

d"'x=p''p'' ■■■p'-diA2---di^ (61) 
and the result of applying evaluated on cr(t). The metric tensor 

5 = -gO e° + e^^ (62) 

is now given in terms of the above as a Taylor series to order A^, and 
generalises the results of [21 UHl EOl EU • 

5 Identities from Gauge Covariance 

Attention is restricted to background gravitational and electromagnetic fields. 
The former is described in terms of a metric tensor g and a metric-compatible 
connection V. The latter is given in terms of a 2-form F. The basic gravita- 
tional variables will be a class of arbitrary local orthonormal 1-form coframes 
{e°} on spacetime related by 50(1,3) transformations. In such frames the 
connection V is represented by the spacetime 1-forms {u;°'^}. The electro- 
magnetic U{1) connection can be represented locally by a spacetime 1-form 
A and F = dA. All other (matter) fields are denoted collectively by ^. It is 
assumed that in a background electromagnetic and Einstein-Cartan gravita- 
tional field described by {e,uj} the matter field interactions can be derived 
from some action functional: 

S[e,uj,A,<i>] = [ A (63) 
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where the 4-form A is constructed from the field variables {e,uj, A,^} and 
their derivatives, so that S is invariant under change of local U{1) and 
Lorentz group transformations ^. Its invariances impose restrictions on the 
variational derivatives that arise in the expression 



Jm Jm 



(64) 



where the vector field X generates a local spacetime diffeomorphism and Cx 
denotes the Lie derivative with respect to X. Since Cx^ = ix dA + dix-A 
and dA = the left hand side of (|64() is zero if X has compact support. 
The 3- forms t^,S^'^, and j are the basic variables in our balance laws and 
may be termed the source currents for gravitation and electromagnetism. In 
the absence of electrically charged matter and a background electromagnetic 
field a recent derivation of the identities that arise from gauge covariance 
can be found in The inclusion of electromagnetic interactions demands 
only a minor modification to these arguments but gives rise to additional 
terms. A sketch of the arguments follows. 

Since the matter field equations are to be satisfied, £ = 0. Using the 
identity 

CxA = ixdA+ dixA (65) 

and recognising that under a C/(l) gauge transformation A A + df for 
some spacetime scalar field / one may write 

CxA = ixF + 6u^i)A (66) 

in terms of the generator 6u{i) of ^ gauge transformation. In a similar 
way with the aid of the structure equations (|36|) . (|37j) the terms Cxe^ and 
^x^^v can be written in terms of the generators of local Lorentz transfor- 
mations and exterior covariant derivatives ^23j : 

/ CxA= [ {T^AixT'^ + T^ADiixen + S^'^ AixR^. (67) 
Jm Jm 



+j A ixF + j A 5u(i)A -T^A 5so{3,i)e^ - 5^,^ A Sso{i,3)^^i^)- (68) 

By gauge invariance of the action Jj^^ A 

I ('^i/(i)+'5so(i,3))A = (69) 

J M 



^In theories with gravitational scalar fields it is possible to extend the local gauge 
symmetries to include Weyl scale covariance and an associated connection 1-form. 
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so the last three terms in (|68|) are zero. If the arbitrary components ix^^ 
of X have compact support one may replace the second term in ()68() by 
Dt^ a ix^'^ and conclude: 

/ (Dr, + T^Afx„r^ + VAix.ii^+jAix.F)W° =0 (70) 
Jm 

where {^a} is the dual frame and {W^} a set of test 0-forms. Similarly 
invariance of the action J^j A under Sso{i,3) alone gives: 

f (I)5/-^(T,,Ae'^-T'^Ae^))>V^ = (71) 

for a set of test 0-forms {W^*}. For smooth tensors on a regular domain of 
M one deduces the local identities: 

Br^ + A ix^ + A ix^ R^^ +jAix^F = (72) 



DS^" = -{r^Ae'' -r'' Ae^). (73) 
In these expressions 

Dr" = dT"+a;% Ar^ 



DS""^ = d S"'^ + A S^f^ + o;^^ A 5"^. 
Finally invariance of the action under Su{i) alone yields: 

/ duii)A= [ jA6uii)A= [ jAdf= [ fdj = 0. (74) 
Ja4 Jm Jm Jm 

For smooth j on a regular domain 

dj = (75) 

since /, with compact support, is arbitrary. 

It should be stressed that (f7^ . (|7n|) . (f7^ are identities in a field theory 
satisfying the conditions used in their derivation. However in a dynamical 
scheme that attempts to approximate the currents ra,5^^,j they provide 
powerful constraints. Thus H75|) is clearly a conservation law (in an arbitrary 
spacetime background), since for any 4 dimensional domain 7W C M with 
dA4 = Tii — T,2 + C and Si, S2 disjoint spacelike hypersurfaces in M, and j 
regular on A4: 

= [ dj= [ j = [ j - [ j (76) 
Jm JdM Jt.1 Jt,2 

if j^j = 0. The integral J is the electric charge content of j in Si. 
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6 Models for the Source Currents 



The balance laws (|72p . ^TS\i . ((75)) above are devoid of predictive content with- 
out further information that relates the source currents to each other and 
the background fields. Such equations of state serve to define the dynamical 
variables in the theory and reduce the number of degrees of freedom sub- 
ject to temporal evolution. The established laws of Newtonian and special 
relativistic field-particle interactions offer valuable guides since they should 
be recovered in appropriate limits. Indeed an effective source model should 
enable one to recover the framework of non-relativistic continuum mechanics 
from the balance laws in a Minkowski spacetime free of curvature and tor- 
sion. More generally it is known how to model a spinless pressureless fluid 
in a pseudo-Riemannian spacetime with curvature and recover the motion 
of dust as geodesies associated with the Levi-Civita connection [HSIISI]- In 
spaces with curvature and torsion analogous equations for spinning hyper- 
fluids can be constructed (35t l36t l37. .38. .39. .40) . 

In J2] Dixon considered the possibility of using the pseudo-Riemannian 
balance laws to construct a multipole analysis of the dynamics of an ex- 
tended relativistic structure. Following work by Pryce he gave cogent 
arguments for a type of subsidiary condition that eliminated peculiar mo- 
tions of electrically neutral mass monopoles in the absence of gravitation. 
The approach advocated here is based on a similar description of a parti- 
cle in terms of the evolution of appropriate multipoles along the timelike 
worldline a. The structure of the particle is therefore given in terms of 
components of the source currents in the Fermi frame Thus we seek a 
consistent dynamical scheme that will determine these components and the 
history a where the balance laws are satisfied to some order in a suitable 
multipole expansion. 

The basic assumption is that the source currents for a classical particle 
are localised on the tubular region U about a. Thus if C U is the 
hypersurface t = constant then the test functions in (f7^ . (|7H|) . (|7^ are 
rapidly decreasing as a function of the geodesic distance A in the 3-ball 
centred on (T{t)'^ . The precise structure of the stress-energy 3- form currents 

defines the nature of "particle" variables used to pass from a field to a 
particle description. The corresponding choice of stress-energy tensor field 
T = •k~^Ta®e°' is analogous to a choice of constitutive relation in continuum 
mechanics relating configuration stresses and powers to strains and work. 
The operator * here denotes the Hodge map^ with respect to the metric g. 
The choice should naturally have the proper Newtonian limit for a spinless 
particle. We shall also ensure that the Matthisson-Papapetrou equations 
for electrically neutral spinning particles emerge in a torsion-free pseudo- 

^To lowest non-trivial order the calculations below are insensitive to the precise fall off 
with A 

^In terms of an orthonormal basis the volume form ★I is chosen as e'' A A A e'^. 
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Riemannian spacetimes 1411 131 . Thus in the absence of self- fields (see later) 
consider: 



r° = e° + xF""^ ^pu* (77) 



j = P*e'. (79) 

in terms of the scalar components P", p onU and constant x- 

Since the background fields are assumed regular on a the components of 
the curvature, torsion and F as well as P", S"^, p can each be expanded as 
a Taylor series of the form (13). The coframe e" and connection forms cj^/j 
are also available as expansions in A so the coefficient of the test function in 
the integrand of each balance can be expressed as a power series in A with a 
multinomial dependence on To extract dynamic information one 

now truncates each series to some order and integrates over St. This may 
be achieved by writing 

= sin 6 cos (j), p^ = sin 6 sin p^ = cos 6 

and noting that if the positive integers ni,n2,n3 are all even then 



{p-r {vT dn = \\ ':jJ:: ^ (80) 



If any integer is odd then the integral is zero. Since 

Ae^ Ae^ Ae^ = X^smOde A dcpA dXA d{ct)+0{X^) 

the lowest order of relevance includes terms up to order A^ . The resulting 
equations yield ordinary differential equations for the dynamical variables 
E"^,P" since the conservation equation dj = requires q = f^^j to be a 
constant. It is convenient to introduce new variables: 

pi _ pi 

pO _ po 

Si — Ti^iik 
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Ei — F^i 

and denote ordinary derivatives with respect to ct by an over-dot. 
In terms of the spin 2-form 

S = e A e'^ 

and the form V = one has iy ★ 5 = -k{V A S). Since V restricts to the 
tangent to a one readily recognises that s = s^e^ = —iv * 5 is dual to 
the Pauli-Lubanski spin vector on a where S = /lAe^ + iy + s and the 
1-forms h = hj^e^ and s are spatial: ix^ h = ixg s = 0. Similarly with 
P = P°Xo + P'' Xk one has on a: 

ipS = (hkPk) e° - P° /i + eijk s'p> (81) 



Then with respect to a non-rotating Fermi frame (17^ = 0), the equations 
(t72|l. ([73|). (f75|) give: 



s = h X A-h |(B X s-hE X h) (82) 



1 y . 

h = A X s - -P + |(B X h-Fs X E) (83) 



po ^ _ A . p + ;^[(h . E) + A • (B X h) + f%k (s X E)fc+ 
+ B, Sj Poi - Poj (B • s) - Ej hk f^ok ] 



(84) 



k 



P ) , =R^f,ko - Pp - (A)fc P° - g(E)fc+ 
+ x[(Bxh), + (A)fe(E.h) + 

+ W f\j - (Ex s)j f + (B X h), f ^ofe - (E • h) f\k 

-f^kn BnSj+Pkj (B-s)]. 

(85) 



17 



7 Involution and Constitutive Relations 



For given external fields E, B, f "^3^ , equations ^ 

are 10 ordinary differential equations involving 13 functions (s,h, A, P"^) 
of t and the constants X; Q- Additional information is required to render 
these equations useful. Since radiation reactions due to self fields have been 
neglected such information should provide a means of determining the ac- 
celeration A of the curve a in terms of other variables, that may themselves 
satisfy a differential system. Since A^ = {Va(^)k, tbe equation for a can 
then be determined in any coordinate system. For consistency the final set 
of equations should reduce to an involutive differential system. So, suppose 
^^(t) denotes collectively the variables (s,h, P") and the background fields 
whose t dependence on a is supposed known. Write (jHSI), ^ 
and the time derivatives of the background field components as 

iA = WAi^,A,x,q). (86) 

Let 

/CB(e,g,x)=0 (87) 

be a set of scalar subsidiary conditions (constitutive relations indexed by 
B) that is appended to the above system. Then the requirement /Cg = 
implies 

^%^H^.(e,A,x,.)k = (88) 

where the relation (|87j) and implicit t dependence of the background has 
been used. If, for any x, Q, (|88() determines A = A(^, q, x) one may insert 
this in (|86|) and solve for given initial conditions. Such a system then 
determines the acceleration A. However ()88|) may only determine some of 
the components of A, or the solution may depend on the nature of the 
external fields or the values of the parameters x^Q- If no suitable choice can 
be found for generic external fields then /C^ = is an unsuitable condition 
for the source model under consideration and a better condition (or model) 
is required. Another possibility is that (|88|) leads to new conditions, say 
Lb'{C,XjQ) = independent of A. In this case L^' = requires: 

^WA{^,A,X,q)\L,K=0 (89) 

These equations are now analysed in the same way as (|H5)) and the procedure 
repeated until A is determined for a consistent X) Q, or the process exhibits 
a manifest inconsistency. 

For models with more parameters and more components in the source 
currents the strategy is similar. Finding a constitutive relation that ren- 
ders a differential-algebraic-system involutive is, in general, non-trivial and 
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further recourse to weak field approximations may be necessary to make 
progress. The above strategy is based on the premiss that radiative or self- 
field contributions to the balance laws are not relevant. Such contributions 
are expected to yield terms containing higher derivatives of A (that may 
be representations of a differential-delay system) and the above involution 
analysis must be generalised accordingly. The inclusion of such terms is a 
challenging and important issue that will not be pursued further here. 

To illustrate this procedure equations (jBHI), ijHll) . will be sup- 
plemented with the relations: 

(r„ - xFj * e^) A ^S'^^ = 0. (90) 

With the source model above these become the Tulczyjew-Dixon conditions: 

Pa S"^ = (91) 

or 

ipS = (92) 

in terms of the spin 2-form. In terms of the particle variables on a, ()92() 
gives: 

hP° = sxP (93) 
and this implies 

P ■ h = (94) 

s • h = 0. (95) 
By contrast the condition 

ivS = (96) 
gives h = . 

An immediate consequence of from (|H^ . (|HH|) is: 



J^afS^^Pj = (s • s - h • h) = (97) 
in any background. 
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8 Special Cases 



It is of interest to analyse (jSU, (|H5|) . (|5H) . in Minkowski spacetime 
devoid of curvature and torsion. First consider the case where x = 0- Then 



s = h X A (98) 



h = Axs-^P (99) 



P° = -P • A (100) 

P = -P°A-gE. (101) 
These immediately imply: 

(s • s - h • h) = h • P (102) 



(-(P°)2 + P-Py= -2gP-E. (103) 

If the spin current is zero, s = and h = (so S^^^ S"^ = 0), then the 
above equations reduce to 

P = (104) 

P° = (105) 

and 

^ = -Q§^ (106) 

without further conditions. Introducing the constant of motion = m the 
covariant equation of motion for the particle follows from the last equation: 

Vaa = (107) 

m 

in terms of the Levi-Civita connection V with a = g~^(a, — ) for any 1 form 
a. One recognises the equation of motion of a spinless particle with mass m 
and electric charge q jl5) . If F = the particle follows a geodesic of this V. 
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Suppose next that T,af3 ^ 0. The early descriptions of spinning point 
particles (see also employed conditions leading to h = or 5 = 0. 

With the equations (jMl), 11001), itTnTT) this yields: 

s = (108) 



P = 2A X s (109) 

P° = (110) 

P = -P°A-gE. (Ill) 
Then 

(s • s) = (112) 

(P„P°y= -2gP -E (113) 

and s = So 7^ and = m ^ are constants of the motion. Differentiating 
((TU^ gives 

2A X So = -mA - gE. (114) 

If g = this implies A • sq = and hence A • sq = 0. Taking the cross 
product of (|114|) (with q = 0) with Sq then gives: 

A = noxA (115) 

where f2o = —or^^o- Thus, in addition to geodesic motion given by A = 0, 
if A (to) 7^ for any Iq, solutions exist where the acceleration exhibits a 
rotation with angular velocity CIq: 

A{t) = n{t) A{0) (116) 

where 7?. is a time dependent 5*0 (3) matrix describing a rotation with 
angular velocity Qq. Such a motion is unnatural for a classical particle 
free of electromagnetic and gravitational interactions. The existence of such 
solutions is one of the reasons that the Tulczyjew-Dixon condition is to be 
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preferred over (|96|) . If g 7^ then the additional terms yield the modified 
equation 

A = fio X A - — E X fio - 4—^0 (E • fto) ( (117) 
m m \ J 

which also exhibits unnatural motions requiring a specification of an initial 
acceleration. 

With ^ instead of ^ equations (jlHI), (jlHI), itTUTl) reduce to 

a system that determines A rather than A and hence fixes the motion of 
the particle in terms of its initial position and velocity. Differentiating the 
condition and using equations (|U5|) . (|^ . (|1U(J|) . (|1U1|) yields a new 
(vector) condition: 

P°P = 2gs X E (118) 

and hence s • P = and E • P = 0. The invariant Pa i-*" is a constant of 
the motion which we take to be the non-zero constant — m^. Differentiating 
H118|) and using the previously obtained equations gives rise to the equation: 

(P°)^ A + (P°)3 qE + 2q {P^f s X E - V ( ^ • s) E x r = (119) 

where r = s x E. For a generic E, (r 7^ 0), one may take r, s,E as a basis 
in which to evaluate A from (|119j) in terms of P'^,s, E,E. Thus 

A = Arr + AsS + AeB (120) 

where 

^. = -^r.(Exs) (121) 



^E = ^A-E-^^-^A-s (122) 



(B-s) - IEP 
A, = -^^^A-E + ^A-s (123) 



and 



A-E=-^-^E^(sxe) (124) 



A.s=-, , (125) 



(p0)2 + 2g (E • h) 
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with 



s X (s X E). (126) 



(PO) 

The above expression for A{P^, s, E, E) can now be inserted into H98() . 1)100(1 
yielding: 

pO = -^E.(Axs) (127) 



2o 

^ A X (s X (s X E)). (128) 



(PO) 

Given initial values for s and solutions of the coupled equations H127() 
and ()128|) can be used to determine the acceleration from (jl2()|) . Since h = 
and s = if E = we set s|^^q = sq constant. Thus the constant of motion 
s-s — h-h= I Sop identifies sq as the classical spin of the particle in the 
absence of electromagnetic interactions. Similarly m may be identified with 
the particle's rest mass in the absence of such interactions. Thus at any 
t = to: 

{P'{to)f =m' + |s X E|2(to) (129) 



Hto)\' = [sop + |(s X (s X E))P(to). (130) 

If the particle moves slowly with velocity v in an inertial Minkowski 
frame, where the electromagnetic field has components £ = 0,B, [IHl, then in 
the Fermi frame of the particle E ~ ^ x ^B, B ~ ^B. Thus the spin coupling to 
the electromagnetic field indicates that it does not have a magnetic moment. 
Such an interaction arises from the inclusion of the terms in the current that 
depend on x- To see this consider x « 1- Then the equations (|98|) . (|99|l . 
(fTUn|) . (fTUH) are replaced by: 

h = -^ + Axs-|(^Exs + hxB) (131) 



s = hxA + |(Bxs + Exh) (132) 
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pu^_p.A + ^(E.(Axs)-^P-E) + 0{x') (133) 



rE-P'JA+x ( B X (A X s) + A (E • h) - i B X P ) +0{x^). (134) 



Consider static fields (E = 0, B = 0). Then after differentiating the 
condition (|92() , taking scalar products of the result with s and P and using 
the above equations along with the relations P • s = and s • h = an ad- 
ditional constraint arises that for P*^ 7^ requires P = at this order. This 
in turn implies (for s 7^ 0) h = and hence = m constant. Furthermore 
now: 

A = -^E (135) 
m 



s = |Bxs. (136) 
But since h = 0, A x s = x (E x s)/2 or 

X = (137) 

m 

Thus, in this small X) stationary field, limit one has (|107)) and 

= \ hf.S'^p (138) 

and one may identify the interaction of electromagnetism with a charged 
particle possessing a magnetic moment related to its spin with a gyromag- 
netic ratio 2 [3^1. Thus the source currents ()142l) . ()143() . (|144() . are appro- 
priate for the description of a spinning particle of mass m, electric charge q 
and gyromagnetic coupling defined by (|137() . 

9 Generalisations 

A relativistic description of a material rod involves a 2-dimensional immersed 
timelike submanifold S of spacetime (the rod worldsheet) with attached ma- 
terial body frames that relate dilation, bending, shear and torsional defor- 
mations to the immersion. As before, to obtain such an immersion from 
an interacting field description requires a reduction, or truncation, of field 
degrees of freedom. This is most naturally accomplished in terms of a field 
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of local adapted ortho-normal spacetime frames generated by the transport 
of local material body frames along spacelike autoparallels of the Lorcntz 
connection V to points in a neighbourhood lA of the worldsheet S. Di- 
mensional reduction of the stress-energy and spin Bianchi identities, and 
entropy imbalance relations on ZY, is performed in new Fermi coordinates 
{(To, (Ti, A,pi,p2} (with P1+P2 = 1) adapted to such body frames. Their 
definition in terms of the exponential maps is a generalisation of that given 
earlier with S replacing the worldline a. In such coordinates the history of 
the rod E is given as the timelike submanifold {Xpi = 0, Xp2 = 0} for some 
range of {(To,cri} and the local coframe on U may be written: 

= da' + (t^ (139) 



= ff^i da' +p'^dX + / (140) 

for i = 0, 1 and a, (3 = 2,3, where f^i, f^i are forms dependent on 
cr*. A, and with 1-forms (p^^ {A = 0,1,2,3) independent of da', dX but 
components dependent on a' , X, p". Similarly in such Fermi coordinates, 
one may write the connection 1-forms in this coframe 

oo^B = A^Bi da' + C^B (141) 

where A'^Bi are forms dependent on a,' X, and C^bi are 1-forms 
independent of da'\ dX but with components dependent on cr*, A, p". All the 
functions appearing in these expressions can be found as before by solving 
the structure equations that define the torsion and curvature of V. 

The source currents ta, S^^ and j are now localised on the tubular 
region U about the rod history S. The parameterisation of the stress-energy 
3- form currents r^, 5*"^^ on S defines the rod variables used to pass from a 
field to a " Cosserat rod " description ^. In the absence of dissipation one 
can generate relativistic Cosserat models interacting with gravitational and 
electromagnetic fields from: 

= P^,* + xF^^ Sbb'* e^' (142) 



= SA^j*e^' (143) 



j = p*e° + J*e\ (144) 

in terms of the scalar components P^i, T^^^j, F^^ , p, J onU and constant 
X- For thermally active media one must supplement these equations with 
phenomenological consitutive relations for the stress-energy, spin, entropy 
and heat forms subject to the constraints imposed by entropy production, 
''strings arise as a special case of shear- free rods 
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10 Conclusions 



With the aid of the structure equations for a metric-compatible connec- 
tion on the bundle of analytic orthonormal frames over a tubular domain of 
spacetime, a method for constructing a local section by Taylor series has been 
described. Explicit computations have been presented in terms of Fermi co- 
ordinates for this domain associated with a (generalised) Fermi- Walker frame 
attached to a timelike curve. This tubular section is used to effect a trunca- 
tion scheme based on balance laws associated with identities derived from 
the invariances of an action integral. The curve is defined as the history of a 
massive particle with spin. A series of auxiliary constitutive relations is used 
to construct models for the stress-energy and electromagnetic currents that 
enter into these balance laws. Taylor series representations are expressed 
in terms of expansions in the arc distance along spacelike autoparallels as- 
sociated with the Fermi- Walker frame. By suitably truncating these series 
one may construct models of spinning matter. In the absence of self- forces 
the determination of a consistent dynamical scheme for the acceleration of 
the particle worldline can be effected by ensuring that the auxiliary condi- 
tions render each truncated set an involutive system of ordinary differential 
equations. This procedure has been demonstrated by analysing the dynam- 
ical equations in Minkowski spacetime where it becomes possible to identify 
the expected interactions of a charged spinning particle with a particular 
magnetic moment in an external electromagnetic field. 

A number of interesting features emerges from this investigation. The 
inclusion of charged matter interactions with gravity yields terms that may 
be interpreted as tidal forces that couple components of the torsion ten- 
sor to the electromagnetic field via the particle spin. Such interactions are 
novel and may offer new experimental probes for the detection of spacetime 
torsion. Another feature is that in the lowest approximation considered 
here, electrically neutral massive spinless particles are predicted to follow 
autoparallels of the Levi-Civita connection. By contrast the precession rate 
and motion of spinning gyroscopes in gravitational fields with torsion dif- 
fer from those in a torsion- free environment [IHl- The differences can be 
estimated from the results given in this paper. 

The approach adopted here offers a number of natural generalisations. 
Fermi charts adapted to spacetime submanifolds with co-dimension less than 
three permit the methods to be applied to derive the relativistic motion 
of Cosserat media in general spacetimes. The associated orthonormal 
coframe and source currents for a relativistic rod have been mentioned in 
IX. 

The above analysis has been conservative. Prompted by puzzles in cos- 
mology many modern theories of gravitation invoke scalar fields that couple 
directly to spacetime geometry and Yang-Mills fields. Some of these ap- 
proaches suggest that at some scale the dilation group may arise as an ad- 
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ditional local symmetry. If this group is included it is possible to generalise 
the balance laws to incorporate fields with Weyl charge in non-Riemannian 
geometries. It is then possible that not all electrically neutral spinless parti- 
cles have the same type of worldline. In geometries with torsion, autoparallel 
worldlines are in general distinct from geodesies of the (torsion-free) Levi- 
Civita connection j44j . Thus the shift in perihelia of celestial objects in 
highly eccentric orbits may also be used as probes in generalised theories of 
gravitation [T5| . 

A discussion of the significance of radiative forces on the motion of 
charged particles in geometries with torsion and an extension of the methods 
to retarded null coordinates associated with an arbitrary timelike curve will 
be presented elsewhere. 
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